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Evidence for an anomalous like-sign dimuon charge asymmetry 
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We measure the charge asymmetry A = (N ++ — N~~)/(N ++ + N~~) of like-sign dimuon events 
in 6.1 fb _1 of pp collisions recorded with the DO detector at a center-of-mass energy ^fs — 1.96 TeV 
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Studies of particle production and decay under the re- 
versal of discrete symmetries (charge, parity and time) 
have yielded considerable insight into the structure of 
theories that describe high energy phenomena. Of par- 
ticular interest is the observation of CP violation, a phe- 
nomenon well established in the K° and B® systems, but 
not in the B° system where the effects of CP-violation 
are expected to be small in the standard model (SM) 1]. 
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A review of the experimental results and of the theo- 
retical framework for describing CP violation in neutral 
mesons decays can be found in Ref. |2|. The violation 
of CP symmetry is a necessary condition for baryogene- 
sis, the process thought to be responsible for the matter- 
antimatter asymmetry of the universe [3|. However, the 
observed level of CP violation in the K° and B% sys- 
tems is not sufficient to accommodate this asymmetry, 
suggesting the presence of additional sources of CP vio- 
lation beyond the SM 0]. 

This Letter and a more detailed Article Q present a 
measurement of the charge asymmetry for like-sign muon 
pairs. The data, corresponding to an integrated luminos- 
ity of 6.1 fb^ 1 , were recorded with the DO detector [6J at 
the Fermilab Tevatron proton-antiproton (pp) collider, 
operating at a center-of-mass energy of 1.96 TeV. The 



DO experiment is well suited to the investigation of the 
small effects of CP violation because the periodic rever- 
sal of the DO solenoid and toroid magnetic field polarities 
results in a cancellation of most detector-related asym- 
metries. In addition, the pp initial state is a CP eigen- 
state, and the high center-of-mass energy provides access 
to mass states beyond the reach of the B-factories run- 
ning at v/i = M(T(45)). 

The like-sign dimuon charge asymmetry A is defined 
as 



A = 



N++ - N- 
N++ + N- 



(1) 



where N ++ and N represent the number of events 
in which the two muons of highest transverse momen- 
tum, satisfying the kinematic selections described below, 
have the same positive or negative charges. After remov- 
ing contributions from background and from remaining 
detector effects, any residual asymmetry is assumed to 
arise solely from the mixing of B q (q = d, s) mesons (via 
B q -R* B® oscillations) that later decay semileptonically. 
This corrected asymmetry A^ is defined as 



Ab _ ^b 
/1„ 



n: 



N? 



n: 



(2) 



At the Fermilab Tevatron collider, b quarks are pro- 
duced mainly in bb pairs. In like-sign dimuon events, 
one muon arises from direct semileptonic decay, e.g., 
b —¥ /J,~X, of a By or B~ meson, and the other muon 
from a B? O B? oscillation followed, in this example, by 
a semileptonic decay of the B q meson, B q — >• B q — >• n~X. 
The main background for these measurements arises from 
events with at least one muon from kaon or pion decays 
or from the sequential decay of b quarks, b — > c — > /i + A. 
The most important background asymmetry arises from 
the fact that K + and K~ mesons interact differently with 
the material of the detector, and thus their decay rates 
into positive and negative muons are not identical. 

The asymmetry A b sl can also be obtained from the mea- 
surement of the charge asymmetry a^ in semileptonic de- 
cays of b quarks to muons of "wrong charge" , i.e., a muon 
of charge opposite to the sign of the charge of the original 
b quark, induced through B q <rr B q oscillations [121 ]: 



b _ 



T(B -> B -> ii+ X) - T(B -► B -> prX) 
T(B -> B -> n+X) + T(B ->• B ->• n~X) 



- A b 



The asymmetry a^ can be measured from the inclusive 
muon charge asymmetry 



(7) 



where N b + 



and N b represent the number of events 



containing two 6-hadrons decaying semileptonically into 
two positive or two negative muons, respectively. Assum- 
ing CPT invariance, each neutral B q meson contributes 
a term to this asymmetry 



A b sl = p d a d sl + p s a: 



sh 



with a q sl 



Ar g 

AM 



tan< 



(3) 

(4) 



where <p q is the CP-violating phase, and AM q and 
AT q are the mass and width differences between the 
eigenstates of the propagation matrices of the neutral 
B° q mesons. The values of f3 d = 0.506 ± 0.043 and 
p s = 0.494 ± 0.043 are taken from previous measure- 
ments 0| . The SM predicts S = 0.0042 ± 0.0014, 



-o.09i±8;811, Ar d 



(26.7l6;5)xl0- 4 ps- 1 ,Ar s 



(9.6 ± 3.9) x 10" 2 ps- 



(-4-8Ii.2) x !0" 4 and 



a s A = (2.06 ± 0.57) x 10" 5 1]. Deviations of the ex- 
pected asymmetries from zero are negligible compared 
to the present experimental sensitivity, and correspond 



to a small value for A h sl [1|: 

4 b ,(SM) = (-2.3 



+0.5N 
OS) 



x 10" 



(5) 



Extensions of the SM containing additional contributions 
to the Feynman "box" diagrams responsible for B® mix- 



ing can result in larger values of cj> q or Ar g [7Hll1| . Mea- 
surements of A^ , cj) q or AT q that differ significantly from 
the expectations of the SM would therefore be indicative 
of the presence of physics beyond the SM. 



where n + and n~ correspond to the number of positive 
and negative muons satisfying the kinematic selections. 
For the asymmetry a, the signal comes from B q mixing, 
followed by the semileptonic decay. In addition to the 
background already considered for A, the direct produc- 
tion of c quark pairs followed by their semileptonic decays 
constitutes an additional source of muons contributing 
to a. 

We define all muons from weak decays of b and c quarks 
as signal, and use the branching fractions and momen- 
tum spectra of particles in the decay chains that produce 
such muons to determine the dilution of the ^4^ asym- 
metry in the observed asymmetry of the signal compo- 
nent. The dilutions, defined as the coefficients which re- 
late the signal asymmetries to A b v are 0.070 ±0.006 and 
0.486 ± 0.032 for the inclusive muon and for the like-sign 
dimuon signal asymmetries, respectively. The difference 
in the dilution coefficients arises because the presence of 
the second muon with the same charge acts as a flavor 
tag. Therefore, the asymmetry A is far more sensitive to 
.A*! than a. 

We measure the asymmetries A and a in the like-sign 
dimuon and inclusive muon data, respectively. These 
have different contributions from background processes 
and from detector asymmetries, which are measured di- 
rectly in data as a function of the muon transverse mo- 
menta, and are used to correct the measured asymme- 
tries. After applying all corrections, the only expected 
source of residual asymmetry in both the inclusive muon 
and dimuon samples is the asymmetry A^. Given the dif- 
ference in sensitivity between A and a and the fact that 



the asymmetry a is dominated by detector effects, we do 
not take a weighted average of the two determinations of 
A^. Instead, we use the measurement of a to constrain 
the background contributions to A, thereby achieving a 
further reduction of the total uncertainty on A^. This 
is possible because the detector effects and their related 
systematic uncertainties largely cancel in the linear com- 
bination of A and a. 

The inclusive muon and like-sign dimuon samples are 
obtained from data collected with single and dimuon trig- 
gers, respectively. Charged particles with transverse mo- 
mentum in the range 1.5 < pr < 25 GeV and with pseu- 
dorapidity \rj\ < 2.2 [13| are considered as muon candi- 
dates. The upper limit on p? is applied to suppress the 
contribution of muons from W and Z boson decays. To 
ensure that the muon candidate passes through the de- 
tector, including all three layers of the muon system, we 
require either px > 4.2 GeV or a longitudinal momentum 
component \p z \ > 6.4 GeV. Muon candidates are selected 
by matching central tracks with a segment reconstructed 
in the muon system and by applying tight quality require- 
ments aimed at reducing false matching and background 
from cosmic rays and beam halo. The transverse impact 
parameter of the muon track relative to the reconstructed 
pp interaction vertex must be smaller than 0.3 cm, with 
the longitudinal distance from the point of closest ap- 
proach to this vertex smaller than 0.5 cm. Strict qual- 
ity requirements are also applied to the tracks and to 
the reconstructed pp interaction vertex. The inclusive 
muon sample contains all muons passing the selection re- 
quirements. If an event contains more than one muon, 
each muon is included in the inclusive muon sample. The 
like-sign dimuon sample contains all events with at least 
two muon candidates with the same charge. These two 
muons are required to have an invariant mass greater 
than 2.8 GeV to minimize the number of events in which 
both muons originate from the same b quark. 

Muons from decays of charged kaons and pions and 
from incomplete absorption of hadrons that penetrate 
the calorimeter and reach the muon detectors ("punch- 
through" ) , as well as false matches of central tracks to 
segments reconstructed in the outer muon detector, are 
considered as detector backgrounds. We use data to mea- 
sure the fraction of each source of background in both the 
dimuon and inclusive muon samples, and the correspond- 
ing asymmetries. Data are also used to determine the 
intrinsic charge-detection asymmetry of the DO detector. 
Since the interaction length of the K + meson is greater 
than that of the K~ meson Q, kaons provide a positive 
contribution to the asymmetries A and a. The asymme- 
tries for other background sources (pions, protons and 
falsely reconstructed tracks) are at least a factor of ten 
smaller. 

The asymmetry for kaon tracks that are eventually 
misidentified as muons (K — > /i tracks) is measured in 
data using K*° -> K^ir? and 4> -» K + K~ decays. For 
both channels we select muon candidates from the entire 
inclusive muon sample and examine mass distributions 
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FIG. 1: The difference in the number of events for the K + w~ 
and K~tv + mass distributions of K*° candidates in the inclu- 
sive muon sample. The solid line represents the result of the 
fit, while the dashed line shows the background contribution. 



separately for events with K + — > /i + and K~ — > fj,~ 
tracks, extracting the sum and the difference in the num- 
ber of K*° or (f> meson decays containing positive or neg- 
ative K — > /i tracks. The distribution of this difference 
as a function of the invariant mass of the K*° candi- 
dates is shown in Fig. [T] The resulting asymmetry is 
corrected using simulations |5J for the fraction of kaons 
(~ 6%) that decay prior to being reconstructed. Simi- 
larly, the asymmetry for pion or proton tracks misidenti- 
fied as muons are measured using samples of K$ —t ir + ir~ 
and A — > p-x~ decays, respectively. 

The fraction of muons from kaons is also determined 
from the K*° —> K + ir~ sample. The fraction of all 
kaons arising from K*° decay is taken from the observed 
K* ± —> A"57T ± decays using the assumption of isospin 
invariance, which is validated in data |2j. The probabil- 
ity of identifying the associated tt^ 1 in the K*° decay is 
taken to be the same as in K* ± , as is confirmed by sim- 
ulation. The fractions of pions and protons associated 
with identified muons, relative to the fraction of muons 
from kaons, are estimated using the decays K$ — > ir + n~ , 
<p —t K + K~, and A — > pir~ decays, and the spectra and 
multiplicities of pions, kaons and protons from simula- 
tion. 

After subtracting the muons originating from kaons, 
pions and protons, we find the fraction of muons in 
the inclusive muon sample from prompt sources consti- 
tuting the signal sample (heavy flavor) to be 0.581 ± 
0.014 (stat) ± 0.039 (syst). The signal fraction arising 
from prompt sources for the like-sign dimuon sample, af- 
ter subtracting the contribution from events where one 
or both muons are background, is 0.665 ± 0.016 (stat) ± 
0.033 (syst). 

Tabic III of Ref. 14] lists all significant contributions 
to the dimuon charge asymmetry caused by detector ef- 
fects. The largest of these is w 3%. The reversal of both 
solenoid and toroid magnet polarities suppresses many 
detector effects, reducing thereby any charge asymmetry 
introduced by track reconstruction and muon identifica- 



tion considerably [14| . The small residual reconstruction 
asymmetry is measured using a sample of J/ip — > M + M~ 
decays reconstructed from two central detector tracks, 
with at least one matching a track segment in the muon 
detector. This measurement is performed as a function 
of the muon pt and indicates a residual detector asym- 
metry of order 10~ 3 . 

The uncorrected asymmetries a and A are obtained by 
counting the number of events of each charge in the in- 
clusive muon and like-sign dimuon samples, respectively. 
There are 1.495 x 10 9 muons in the inclusive muon sam- 
ple and 3.731 x 10 6 events in the like-sign dimuon sample. 
The uncorrected asymmetries are 



a = +0.00955 ± 0.00003 (stat), 
A = +0.00564 ±0.00053 (stat). 



(8) 
(9) 



After correcting for background and for the dilutions of 
the j4*j asymmetry in the observed asymmetries of the 
signal component, we obtain 



TABLE I: Sources of uncertainty on A b sl in Eqs. JTDJ, (fill) , 
and (|13p . The first eight rows correspond to statistical uncer- 
tainties and the next three rows to systematic uncertainties. 



Source 



s(A b sl )m 8(A b sl )m s(A b sl )m 



Total statistical 

Total systematic 

Total 



0.01118 
0.02140 
0.02415 



0.00266 
0.00305 
0.00405 



A or a (stat) 0.00066 0.00159 0.00179 

K fraction (stat) 0.00222 0.00123 0.00140 

vr fraction 0.00234 0.00038 0.00010 

p fraction 0.00301 0.00044 0.00011 

K asymmetry 0.00410 0.00076 0.00061 

7T asymmetry 0.00699 0.00086 0.00035 

p asymmetry 0.00478 0.00054 0.00001 

detector asymmetry 0.00405 0.00105 0.00077 

K fraction (syst) 0.02137 0.00300 0.00128 

ir, K, p multiplicity 0.00098 0.00025 0.00018 

At dilution 0.00080 0.00046 0.00068 



0.00251 
0.00146 
0.00290 



A b sl = +0.0094 + 0.0112 (stat) ± 0.0214 (syst) 



(10) 



from the inclusive muon sample and 

-0.00736 ± 0.00266 (stat) ± 0.00305 (syst) (11) 



A b 

Ail 



from the like-sign dimuon sample. The uncertainties 
on the measurement obtained from the like-sign dimuon 
sample are much smaller, as expected from the factor of 
seven difference in the dilution coefficients that relate the 
signal asymmetries to A^ in the two samples. Since the 
same background processes contribute to the uncorrected 
asymmetries a and A, their uncertainties are strongly cor- 
related. We take advantage of this correlation to obtain 
a single optimized value of A b ^ with higher precision, by 
using a linear combination of the uncorrected asymme- 
tries 



A' = A 



aa. 



(12) 



We scan the coefficient a in order to minimize the total 
uncertainty on the value of A^, which occurs when a = 
0.959. The corresponding final result for the asymmetry 
A b sl is: 



A b 

A sl 



-0.00957 + 0.00251 (stat) ± 0.00146 (syst). (13) 



It differs by 3.2 standard deviations from the SM predic- 
tion for A b t of Eq. J5]). The contributions to the total 
uncertainty of A b t in Eqs. ([TO]). (fTT|) and (JTHJ are listed 
in Table HI and the result in Eq. (IT31 is dominated by 
statistical uncertainties. 

Several consistency checks are performed by dividing 
the data into smaller samples using additional selections 
based on data taking periods, muon and track quality 
requirements, and changing the requirements on impact 
parameter, transverse momentum, polar angle and ra- 
pidity of the muons. The resulting variations of A b \ are 
statistically consistent with the result of Eq. Q21 even 
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FIG. 2: (Color online) The observed and expected like-sign 
dimuon charge asymmetry A as a function of dimuon invari- 
ant mass. The expected asymmetry is shown for A al =0 and 
^,=-0.00957. 



if the individual values of the uncorrected asymmetries 
A and a vary widely between the different samples due 
to changes in the background contributions. Both the 
size and the dependence on the muon momentum of the 
asymmetry a, which is dominated by background, are 
reproduced correctly through measurements of the back- 
ground fractions and asymmetries. Similarly, the depen- 
dence of the like-sign dimuon asymmetry on the dimuon 
invariant mass observed in data, as shown in Fig. [2] is 
reproduced by expectations when A^ is fixed to its mea- 
sured value, while there are significant discrepancies if 
A b { = is assumed. 

The measured value of A^ places a constraint on the 
charge asymmetries in semileptonic decays of B^ and B® 
mesons and on the CP-violating phases of the J3° and 
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FIG. 3: (Color online) Comparison of A^ in data with 
the SM prediction for a^ and aj. Also shown are 
other measurements of 4 = -0.0047 ± 0.0046 [H^ and 
o|i = -0.0017 ±0.0091 [if]. The bands represent the ±1 
standard deviation uncertainties on each measurement. 



B® mixing matrices, as given by Eqs. [3]and[4] Figure [3] 
presents our measurement of A^ in the a^-a^ plane, 
together with direct measurements of a^ from the B- 
factories |15l - [l7| and of our independent measurement of 
a|j in B® — > D s /iX decays [lq . Additional comparisons 



and combinations of these results with previous measure- 
ments sensitive to the same physics effect are given in 
Rcf. [|. 

In conclusion, we have measured the like-sign dimuon 
charge asymmetry A^ of semileptonic &-hadron decays: 

A b sl = -0.00957 ± 0.00251 (stat) ± 0.00146 (syst). (14) 

This measurement is obtained from data corresponding 
to 6.1 fb _1 of integrated luminosity collected with the DO 
detector at the Fermilab Tevatron collider. It is consis- 
tent with our previous measurement |14| obtained with 
1 fb _1 and supersedes it. The asymmetry disagrees with 
the prediction of the SM by 3.2 standard deviations. This 
is the first evidence for anomalous CP violation in the 
mixing of neutral B mesons. 
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